Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

30J能量驱动类镍银X光激光的理论研究

 在1.2×10¹³W·cm^－2低功率密度下，对基频激光预主短脉冲驱动类镍银X光激光机理进行了数值模拟和理论分析。证实了在靶长23 mm范围内X光激光都能获得有效放大，取得了和实验相符合的结果。考虑了单柱面镜线聚焦沿靶长度方向功率密度的非均匀性对X光激光放大的影响，采用弯曲靶能有效克服折射以及单柱面镜线聚焦功率密度非均匀带来的不利影响。理论模拟给出的类镍银X光激光的出光下限泵浦功率密度也与实验符合得很好。理论模拟还表明，采用1%左右的预脉冲强度并对预主脉冲时间间隔进行优化，X光激光的输出能量和能量转换效率将获得大幅度提高。     

晶格失配对异质外延超薄膜生长中成核特性的影响

利用动力学蒙特卡罗方法模拟了异质外延超薄膜生长中的成核过程.研究了薄膜与衬底的晶格失配对超薄膜生长中成核密度、平均核尺寸、标度关系及生长模式的影响.结果发现产生压（张）应变的晶格负（正）失配使生长过程更早（迟）从成核区进入过渡区，失配越大，这一效应越明显.在相同的沉积条件下，负失配导致超薄膜形成较低的成核密度与较大的平均核尺寸，而正失配则相反.成核密度满足标度关系Ns≈(F/D)χ，随着失配度从-0.04增加到0.02，标度系数χ从0.37逐渐减小到0.33，对应超薄膜生长过程从包含二聚体扩散模式转变到无 关键词： 薄膜生长 成核 晶格失配 蒙特卡罗模拟     

Cl/HN3/I2全气相化学激光体系的化学动力学模拟计算

 对Cl/HN₃/I₂产生NCl(a)/I激光的过程进行了化学动力学计算，主要考察了Cl,HN₃和I₂的初始粒子数密度及其配比对小信号增益系数的影响。结果发现,当温度为400K, 初始Cl粒子数密度为1×10¹⁵，1×10¹⁶和1×10¹⁷cm^-3时,小信号增益系数分别达到1.6×10^-4,1.1×10^-3和1.1×10^-2cm^-1，获得最佳小信号增益系数的HN₃和I₂的初始粒子数密度分别为初始Cl粒子数密度的1～2倍和2%～4%。同时，对Cl，HN₃和I₂配比对小信号增益系数和增益持续时间的影响进行了讨论。     

19.6nm波长类氖锗X光激光光源理论模拟

 波长19.6nm的类氖锗X光激光适合作为诊断激光等离子体界面不稳定性的光源。用经过实验检验的系列程序对预-主短脉冲驱动类氖锗进行了系统的优化设计和理论分析。采用2%~3%的预脉冲强度，6~8ns的预-主脉冲时间间隔，在4×10¹³W/cm²功率密度驱动下, 波长19.6nm增益区的宽度可以超过60μm，增益区的维持时间可以达到90ps。对于16mm长的平板靶，增益系数可达11.8/cm；弯曲靶增益系数可达13.3/cm；单靶小增益长度积可达21.3，单靶就可以获得饱和增益。采用双靶对接，其小讯号增益可达38.4，可以获得深度饱和增益，能满足应用演示所需的X光激光光源。     

超短脉冲超强激光与固体靶相互作用中高能离子的产生

 用2D3V PIC粒子模拟方法得到了超短脉冲超强激光与固体靶相互作用中高能离子产生的图像，并对其机理进行了研究。在靶前后表面都观察到了高能离子的产生，并诊断了离子能谱。模拟结果表明，在靶前表面所产生的高能离子，角分布较大，在向靶内输运过程中会损失能量；在靶后表面产生的高能离子，定向性很好，能获得很高的能量。模拟得到的离子能量和实验观测结果在量级上相符。     

Studies on the Adsorption of Asymmetrical Triblock Copolymers by Scheutjens-Fleer Theory and Monte Carlo Simulation

The adsorption of asymmetrical triblock copolymers from a non-selective solvent on solid surface has been studied by using Scheutjens-Fleer mean-field theory and Monte Carlo simulation method on lattice model. The main aim of this paper is to provide detailed computer simulation data, taking A8-kB20Ak as a key example, to study the influence of the structure of copolymer on adsorption behavior and make a comparison between MC and SF results. The simulated results show that the size distribution of various configurations and density-profile are dependent on molecular structure and adsorption energy. The molecular structure will lead to diversity of adsorption behavior. This discrepancy between different structures would be enlarged for the surface coverage and adsorption amount with increasing of the adsorption energy. The surface coverage and the adsorption amount as well as the bound fraction will become larger as symmetry of the molecular structure becomes gradually worse. The adsorption layer becomes thicker with increasing of symmetry of the molecule when adsorption energy is smaller but it becomes thinner when adsorption energy is higher. It is shown that SF theory can reproduce the adsorption behavior of asymmetrical triblock copolymers. However, systematic discrepancy between the theory and simulation still exists.The approximations inherited in the mean-filed theory such as random mixing and the allowance of direct back folding may be responsible for those deviations.     

高速色散补偿系统中喇曼放大器性能研究

利用喇曼放大器的非线性偏微分耦合方程组,对带喇曼放大器的色散补偿系统进行了仿真,并与传统的仿真方法进行了比较.结果表明:在色散完全补偿系统中,当传输速率较高时,忽略色散后的传统仿真方法的结果,与采用色散补偿技术后使整个系统色散为零的仿真结果有一定偏差,此时不能忽略色散的作用.还对前补偿系统和后补偿系统进行了研究,利用Q值性能判别法,分别得到了信号光脉冲的最佳占空比和最佳脉冲阶数.